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Abstract 

We present an approach which enables to identify phase synchronization in coupled chaotic 
oscillators without having to explicitly measure the phase. We show that if one defines a typical 
event in one oscillator and then observes another one whenever this event occurs, these observations 
give rise to a localized set. Our result provides a general and easy way to identify PS, which can 
also be used to oscillators that possess multiple time scales. We illustrate our approach in networks 
of chemically coupled neurons. We show that clusters of phase synchronous neurons may emerge 
before the onset of phase synchronization in the whole network, producing a suitable environment 
for information exchanging. Furthermore, we show the relation between the localized sets and the 
amount of information that coupled chaotic oscillator can exchange. 
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I. INTRODUCTION 



The emergency of collective behavior among coupled oscillators is a rather common phe- 
nomenon. In nature, one typically finds interacting chaotic oscillators which through the 
coupling scheme form small and large networks. Surprisingly, even though chaotic systems 
possess an exponential divergency of nearby trajectories, they can synchronize due to the 
coupling, still preserving the chaotic behavior 
have been found in a variety of fields as ecology 



lasers 



3[]. Indeed, synchronization phenomena 
4(, neuroscience {5, 6, 7], economy [gj], and 



In the last years some types of synchronization have been reported 12|. A rather interest- 
ing kind is a weak synchronization, namely phase synchronization (PS), that does not reveal 
itself directly from the trajectory, but as a boundedness of phase difference between the 
interacting oscillators. In such a synchronization the trajectories can be uncorrelated, and 
therefore, the oscillators present some independence of the amplitudes, but still preserving 
the collective behavior. 

This phenomenon can arise from a very small coupling strength 13| . It has been reported 



that it mediates processes of information transmission and collective behavior in neural and 
active networks jl4J ]. and communication processes in the human brain HQ. Its presence 



has been found in a variety o 
in electrochemical oscillators 



experimental systems, such as in electronic circuits 



17, 



191 ] . plasma physics [20j], and climatology [21] . 



181], 



In order to state the existence of PS, one has to introduce a phase 4>(t) for the chaotic 
oscillator, what is not straightforward. Even though the phase is expected to exist to a 



general attractor, due to the existence of the zero Lyapunov exponent [12j, its explicit 
calculation may be impossible. Actually, even for the simple case of coherent attractors, it 
has been shown that phases can be defined in different ways, each one being chosen according 
to the particular case studied. However, all of them agree for sufficiently coherent attractors 



221. 



In spite of the large interest in this field, there is still no general, systematic, and easy 
way to detect the existence of this phenomenon, mainly, due to the fact that the phase is 
rather difficult (often unknown) to calculate. The calculation becomes even harder if the 
oscillators are non-coherent, e.g. the funnel oscillator [121 ] . Therefore, in order to present a 
general approach to detect PS, with practical applications, we must overcome the need of a 
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phase. 

In many cases the phase can be estimated via the Hilbert transformation or a wavelet 
decomposition [12]. Supposing that it is possible to get a phase, the approach developed in 



Ref. 



23] gives rather good results. It is grounded on the idea of conditional observations of 



the oscillators. Whenever the phase of one of the oscillators is increased by 2ir, we measure 
the phase of the other oscillators. The main idea is that if one has PS, the distribution of 
these conditional observation in the phase presents a sharp peak, and therefore PS can be 
detected. 

There are a few approaches that try to overcome the difficulties of not having a general 
phase. For periodically driven oscillators, there is an interesting approach, very useful and 
easy to implement that overcomes the need of a phase, the stroboscopic map technique. It 
consists in sampling the chaotic trajectory at times nTo, where n is an i nteg er and To is the 



period of the driver. The stroboscopic map was used to detect PS [12j, [18|, |2Cj • The basic 
idea is that if the stroboscopic map is localized in the attractor, PS is present. Actually, 
the stroboscopic map is a particular case of the approach of Ref. {23J]. Indeed, since the 
driver is periodic, the observation of the trajectory of the chaotic oscillators at times nTo 
is equivalent to observe the oscillators at every increasing of 2ir in the phase of the driver. 
Furthermore, if the chaotic oscillator presents a sharp conditional distribution, this means 
that the stroboscopic map is localized. The advantage of such an approach is that it does 
not require the introduction of a phase neither in the periodic oscillator nor in the chaotic 
one. 

In the case of two or more coupled chaotic oscillators, namely Sj and the stroboscopic 
map techniques can be no longer applied. However, if the oscillators are coherent and have 



a proper rotation, a generalization of the stroboscopic map has been recently developed [24]. 
Instead of observing the oscillators at fixed time intervals, multiples of the period, one can 
define a Poincare section in Sj and then observe every time the trajectory of crosses 
the Poincare section. If the oscillators are in PS, these observations give place to a localized 
set. 

Another approach that is relevant to the present problem is the one developed in Ref. 



25] . This approach consists of defining a point x.j(t) G and a small neighborhood of 



this point composed by points Xj(tj) G Xj, where i = 1, . . . , N, with N being the number of 
points within the defined neighborhood. Then, one observes the oscillator at the times ti, 
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which gives place to the points Xfc(tj) e E^. Again, the idea is that if the oscillators present 
synchronization, the cloud of points Xfc(t,) occupies an area much smaller than the attractor 
area. Further, estimators have been introduced to quantify the amount of synchronization 
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251 ] . We show that all these 



25]- 

Even though the intuition says that localized set implies the presence of synchronization, 
there is a lack of theoretical analysis showing such a result for a general oscillator. Moreover, 
as far as we know, there are no results that guarantee that such an approach works for 
multiple time-scale oscillators. In addition, it is not clear what kind of points (events) could 
be chosen, and finally, how one should proceed in the case that the small neighborhood of 
the point Xj(£) e Ej has infinitely many neighbor points. 

In this work, we extent the ideas of Ref. [l2, 18, 20 . 
approaches can be put in the framework of localized sets. Our results demonstrate that for 
general coupled oscillators £ 3 - and if one defines a typical event in Ej and then observes 
the oscillator E& whenever this event occurs, these observations give rise to a localized set 
in the accessible phase space if PS exists. These results can be applied to oscillators that 
possess multiple time-scales as well as in neural networks. As an application, we analyze the 
onset of PS in neural networks. We show that in general neural networks one should expect 
to find clusters of phase synchronized neurons that can be used to transmit information in 
a multiplexing and multichannel way. Finally, we relate the localized sets from our theory 
to the information exchange between the coupled chaotic oscillators. 

The paper is organized as follows: In Sec. [IT] we define the dynamical systems we are 
working on. In Sec. IHII we give a result that enables the identification of PS without having 
to measure the phase. We illustrate these findings with two coupled Rossler oscillators in 
Sec. IIVI For oscillators possessing multiple time-scales our main results are discussed in 
Sec. [Vj and then illustrated in Sec. |VT]for bursting neurons coupled via inhibitory synapses. 
Our results are also applied to neural networks of excitatory neurons in Sec. IVHI We briefly 
discuss how to apply these ideas into high dimension oscillators and experimental data series 
in Sec. IVIIII Finally, we analyze the relation between the localized sets and the transmission 
of information in chaotic oscillators in Sec. HXl Moreover, in Appendix |A] we prove the main 
theorem of Sec. II III about the localization of sets in PS. 
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II. BASIC SET UP 



We consider N oscillators given by first order coupled differential equations: 

N 

±i = Fi(Xi) + CijUjiXj-X;) (1) 

3=1 

where, x$ G R™, and Fj : W H — > W H , Hj is the output vector function, and C„ is the 
coupling strength between j and «. Note that could also depend on the coordinates and 
on time. From now on, we shall label the coupled oscillator Xj by subsystem E;. Next, we 
assume that each Ej has a stable attractor, i.e. an inflowing region of the phase space where 
the solution of Ej lies. Further, we assume that the subsystem Ej admits a phase 4>j(t). 
Therefore, the condition for PS between the oscillators Ej and E& can be written as: 

\rmf>j(t) -ruf> k (t)\ < c, (2) 

where n and m are integers, and the inequality must hold for all times, with c being a finite 
number. For a sake of simplicity, we consider the case where n — m — 1, in other words 
1 : 1 PS. Herein, we suppose that a frequency Qj can be defined in each subsystem Ej, such 
that: 

0j = fij(Xx, • • • ,Xjv,t), (3) 

where f2j is a continuous function bounded away from zero. Furthermore there is a number 
M such that f2j < M. This phase is an abstract phase in the sense that it is well defined, 
but we are not able to write the function Qj for a general oscillator. We also consider the 
frequencies 0j not to be too different, such that, in general, through the coupling PS can be 
achieved. 

III. LOCALIZED SETS IN PS STATES 

In this section we present our main result. The basic idea consists in the following: Given 
two subsystems E fc and Ej, we observe E fc whenever an event in the oscillator Ej happens. 
As a consequence of these conditional observations, we get a set V k . Depending on the 
properties of this set one can state whether there is PS. 

The conditional observations could be given by a Poincare section, if it is possible to 
define a Poincare Section with the property that the trajectory crosses it once per cycle in a 
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given direction. We wish to point out that in this case, one is able to have more information 
about the dynamics and the phase synchronization phenomenon. As an example, one can 
introduce a phase, and estimate the average frequency of the oscillators. However, these 



techniques based on the Poincare section 
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24| cannot be applied to attractors without a 



proper rotation, where such a section cannot be well defined. 

Our main result overcomes the need of a Poincare Section. We show that one can use any 
typical event to detect PS. Such events may be the crossing of the trajectory with a small 
piece of a Poincare section (when it is possible to defined such a section), the crossing of 
the trajectory with an arbitrary small segment, the entrance of the trajectory in an e-ball, 
and so on. The only constraint is that the event must be typical (we shall clarify what we 
mean by typical, later on ) and the region where the event is defined must have a positive 
measure. Let (t\j)i^n be the time at which the ith event in the subsystem Sfcj happens. 
Then, we construct the set: 

V k = U ieN F%°), (4) 

where is the initial point within the attractor of £/-. Next, we define what we understand 
by localized set. 

Definition 1 Let Dj be a subset of $ j . The set Dj is localized in if there is a cross 
section and a neighborhood Aj of tyj such that T>j U Aj = 

An illustration of the definition is given in Fig. [TJ 




FIG. 1: (Color online) An illustration of the Def. [TJ The set T>j does not intersect the neighborhood 
Aj, therefore, T>j is a localized set of Qj. 

Under the assumptions of Sec. HIl the following result connects the existence of phase 
synchronization with the localization of sets the T>: 
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Theorem 1 Given a typical event, with positive measure, in the oscillator S 3 -, generating 
the times (t*)j 6 N- The observation ofHk at (ij)ieN generates a localized setT>k if there is PS. 



This result constitutes a direct generalization of approaches of Refs. [2J, |25(. As a 
consequence, this result shed a light into the problem of PS detection, which turned out to 
be a rather difficult task, depending on the system faced. Therefore, PS can be detected 
in real-time experiments and in data analysis by verifying whether the sets T> are localized, 
without needing any further calculations. 



A. Connection between V and Unstable Periodic Orbits 

In this section we investigate the mechanism for the non localization of the sets T>. We 
let the event definition be an entrance in an e-ball in both subsystems, with e being the 
radius. When e is small enough, we can demonstrate that PS leads to the locking of all 
unstable periodic orbits (UPO) between the subsystems. 

Proposition 1 If the set V k is localized, then all UPOs between E fc and Ej are locked. 

Proof: We demonstrate this result by absurd. Let us assume that there is PS; as 
a consequence the set V k is localized. Suppose that there is an UPO, regarded as Xj 
in Ej, and another UPO, regarded as Xj. in E&, and that they are not locked (there is 
no rational number that relates both frequencies). So, there is a mismatch between the 
frequencies of the two UPOs. Given an £j-ball around y° (resp. y°), where y° G Xj (resp. 
y° G Xk), any point x 3 - distant Sj from y°, where 6j <C £j, follows d[F|(xj), F*-(y?)] < Ej, 
for any t < t « £n(ej/5j)/X max , where A max is the largest eigenvalue associated with the 
orbit Xj, and d[-,-] is a metric. An initial condition inside the (^-ball is governed by the 
UPO Xj till a time t < t, see Fig. [2] for an illustration. Next, we construct the set T> k 
by sampling the trajectory of E^ whenever the trajectory Ej enters in the £-ball, which is 
equivalent to observe E^ every period of the UPO Xj. There is an one-to-one correspondence 
(isomorphism) between the dynamics of the conditional observations and the dynamics of 
the irrational rotation in the unitary circle, R a '■ S 1 — > S 1 , R a = e a%/ ^z, where a is the 
frequency mismatch between the two UPOs, here given by: 

a = inf a ,b{axu Xj -bxijj Xk }, (5) 
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radius f'(8j) « e 




FIG. 2: Illustration of the dynamics near a UPO. 

where 0Jx k is the angular frequency of Xj t k- This means that the points of T>k will be 
dense around the UPO Xk, and therefore, the set T>k is not localized; there is no PS, what 
contradict our assumption. Indeed, since Au > 0, it is impossible to bound the phase 
difference between Sj and by a finite number. Thus, in order to have localized sets T>, 
all UPOs must be locked. □ 

This shows that the mechanism for the non-localization of the sets T> will be the existence 
of unlocked UPOs between Ej and Similar results have been pursued for periodically 
driven oscillators, Right at the desynchronization some UPOs become unlocked and 

the stroboscopic map becomes non-localized, and some phase slips happen, generating an 
intermittent behavior. The duration of the phase slips are related to the number of unlocked 
UPOs. Of course, in this regime the set T> is a non-localized set. However, if one looks for 
finite time intervals the set T> may be apparently localized. 



IV. COUPLED ROSSLER OSCILLATORS 

We first illustrate this result for two coupled Rossler oscillators, given by: 

*i,2 = -"1,22/1,2- Zifi + e{x 2> i- x 1>2 ), 

2/1,2 = "1,2^1,2 + 0.152/1,2, (6) 

z 1)2 = 0.2 + z li2 (xi i2 - 10), 

with «! = 1, and a 2 = ot\ + Sa 2 . In such a coherent oscillator, we can simply define a phase 
tancfti = Vi/xi, where i — 1,2, which provides an explicity equation for it. Indeed, taking 
the derivative with respect to time: 

d . , . • d yi ,_. 
g-tan^) x 0, = --, (7) 



which can be written as sec 2 (4>i) x ^ = (2/^ — yix^/xf, which provides: 

Mt)= f "'"'I' T'"- (8) 
Jo x i + vt 

noting that sec 2 (pi = (xf + yf)/x 2 . In a more compact notation, we consider x* = jfe), 
then Eq. (jSJ) can be written as 

Ut) = f ^j^dt, (9) 

JO l X i| 

where A represents the vectorial product. Equation can be used to calculate the phase 
of the oscillators Sj, and there is PS if A<p = <p 2 — 4>i remains bounded as t —>■ oo. 

In order to apply our results we may define an event occurrence in both oscillators. We 
define the event in oscillator Si to be the trajectory crossing with the segment: 

Si = {xi, j/i, zi e R\x 1 < -13, ?/i = 0, and y\ > 0}, (10) 

the crossings generate the times (t\)i € ^. The event in the oscillator S 2 happens whenever 
its trajectory crosses the segment: 

^2 = {x 2 , 2/2, z 2 E R\x 2 > 5, 2/2 = 10, and y 2 < 0}, (11) 

the crossings generates the times (t|)«eN- Then, the set P 2 ,i is constructed by observing the 
oscillators S 2 ,i at times (t\ 2 )ieN 

For e = 0.001 and Aa 2 = 0.001, the set V>\ spreads over the attractor of Si [ Fig [3] 
(a)], and T> 2 spreads over the attractor of S 2 [ Fig. G9^b)]. Therefore, there is no PS, i.e. 
the phase difference A<j) diverges [Fig [3] (e)]. Indeed, a calculation of the frequencies shows 
that {(pi) = 1.03479 and (0 2 ) = 1.03508. As we increase the coupling, PS appears. In 
particular, for e = 0.011 and Aa = 0.001, the sets T>i and D 2 are localized [Figs [3] (c) and 
(d), respectively]. Hence, the phase difference is bounded [Fig. [3(f)] • The average frequency 
is (0i) = (0 2 ) = 1.03522. 

A. Estimating the synchronization level 

Our main goal is to state the existence of PS, however, we can also estimate the syn- 
chronization level between Sj and S& by means of the localized sets. This can be done by 
introducing an estimator Hj^. One way to estimate the amount of synchrony is to define: 

vol of V . 
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FIG. 3: (Color online) PS onset in two coupled Rossler oscillators. In (a,c) we depict the attractor 
of the oscillator Ei and in (b,d) the attractor of £2 hi light gray, the sets T> are depicted in black. 
The bars on Figs, (a) and (c) represent the segment Si, while in Figs, (b) and (d) the segment 
52- In (a) and (b) the sets T>\ and T>2 spread over the attractor of the oscillator Si and £2, 
respectively; and there is no PS, the phase difference diverges (e). The parameters are e = 0.001 
and Aa = 0.001. In (c) and (d) the sets T>\ and T>2, respectively, are localized and there is PS; the 
phase difference is bounded (f). The parameters are e = 0.011 and Aa = 0.001. 



where vol denotes the volume 26]. If there is no PS, the set T>j spreads over the attractor 
of Ej, see Fig. [3](a,b), then, = 1. As the oscillators undergo a transition to PS, Hj k 



27|]. 



becomes smaller than 1. The lower Hjf. is the stronger the synchronization level is 

For attractors with the same topology as the Rossler oscillator, Hjk can be easily calcu- 
lated. Instead of computing the volume, we calculate the area occupied by the attractor in 
the plane (x, y). The area Aj of the attractor of Ej can be roughly estimated by the area of 
the disk with radii r m and tm, see Fig. HJ Thus, Aj = vr(r| f — r^J. On the other hand, the 
set Ej is confined into an angle £ [Fig. H|. Therefore, the area of the set T>j can be estimated 
as i{r 2 M — r^J/2. Thus, the estimator can be written as: 

= i (13) 

We have used Eq. (1131) to estimate the amount of synchronization between the two 



10 



FIG. 4: (Color online) Illustration of a localized set in a Rossler like attractor. The attractor can 
be approximated by a disk with major radius vm and minor radius r m . The sets T> are confined 
within an angle £. 

coupled Rossler of Eq. ([6]). We fix e = 0.001 and vary the mismatch parameter 5a within 
the interval [—0.002, 0.002]. For \5a\ m 0.0009 the coupled Rosslers phase synchronize, which 
means that the set T>j is localized. Therefore, Hjk < 1. The smaller the value of \8oc\ is 
the more localized the set T>j becomes, meaning that the oscillators are more synchronized, 
leading Hj k to low values. At \Sa\ = the two coupled oscillators present their strongest 
synchronization with Hjk ~ 0.22. The results are depicted in Fig. [51 




8a 



FIG. 5: Hjk is depicted for the two coupled Rosslers, Eq. ([6]), with e = 0.001. The estimator Hjk 
is computed by means of Eq. (|13fl . whenever Hij = 1 there is no PS. At \6a\ ~ 0.0009 the coupled 
Rosslers undergo a transition to PS, and therefore, Hjk < 1, which shows the presence of PS. 
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V. OSCILLATORS WITH MULTIPLE TIME-SCALES 



In oscillators with only one time-scale, i.e. one typical period, a typical event means an 
event possible to realize, thus with positive measure. In oscillators with multiple time-scales, 
i.e. oscillators that possess more than one typical period (an oscillator with a fast and slow 
variables), a typical event means an event that takes into account all time-scales. Conversely, 
an atypical event is the one that takes into account just a few time-scales, e.g. only one. In 
such an oscillator with multiple time-scales, one may have synchronization only in one time 
scale, while the others may be asynchronous. If the event definition excludes completely 
the dynamics of the synchronized time-scale this event is atypical and one does not observe 
localized sets through it. In order to clarify these ideas, we consider two instructive examples. 

A. Dynamics on a Torus 

Let us consider a quasi-periodic motion on a torus T 2 with two independent frequencies 
u) and a, i.e. nuj — ma 7^ Vn, m G Z. The dynamics on the torus E& : T 2 — > T 2 
can be characterized by the angular variables and the flow takes the form = (w, v) = 
(out + u>o,ctt + «o). Furthermore, we consider another oscillator on a quasi-periodic torus 
with two independent frequencies u and (3, the flow Ej : T 2 — > T 2 , in angular variables, 
takes the form Qj = (g, h) = {ut + u , (3t + /3 ). Therefore, under this construction one sees 
PS in only one time-scale, since a, (3 are independent. 

If we consider the event in the oscillator E& to be the increasing of 2tc on the variable 
u, conversely the crossing in the section p, it generates the times t\ = 2n x i/uo. The 
observation of Ej at these times generates a localized set Vj, which will lay on S 1 , a subset 
of T 2 = S 1 x S 1 , and will never occupy the full space. On the other hand, if we consider the 
event in the oscillator Ej to be the increasing of ir on the variable h, conversely the crossing 
with the section p, the set T>k will not be localized, since a, (3 and uj are independent. 

Therefore, one must define an event that captures the dynamics of the synchronized 
time scale. In the pictorial example of Fig. [6] any other piece of section that is a linear 
combination of p and p provides typical events. 
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FIG. 6: (Color online) Illustration of two possible sections on the torus T 2 . In (a) the section 
/j, takes into account only the dynamics of u>, the synchronized scale. In (b) on the other hand, 
the dynamics of the synchronized time scale is ruled out on the section p. Therefore, using this 
particular section one cannot observe localized sets T>, since the synchronized time-scale is not 
taken into account. 

B. Spiking/Bursting Dynamics 

An interesting situation is when the time scales present a relationship, which is the case 
for spiking/bursting oscillators. Consider two spiking/bursting neurons A/} and A/jt- They 
have distinct time-scales, the bursting scale, with low frequencies, and the spiking scale, 
with high frequencies. 



The spiking scale consists of the action potentials [28( which occur due to the exchange of 
ions like K + of the external media with the neuron. On the other hand, the neuron may also 
exchange slow current like Ca +2 which inhibits the occurrence of spikes generating the bursts. 
An event denned by the occurrence of a burst defines simultaneously the beginning and the 
ending of a spike train. Therefore, even though spikes and bursts may have independent 
frequencies, the burst occurrence is also determined by the occurrence of the first and last 
spike within the burst. 

It has been reported that it is possible to have PS in the bursting scale while the spiking 
scale is not synchronized 6|. Therefore, in order to analyze the existence of synchronization 
between the neurons, by means of standard techniques, the spiking and bursting scales must 
be separately analyzed. Our method detects PS independently on the time-scale that the 
event is defined; if one time-scale is synchronous one finds localized D sets. In order to 
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Q- We define a 



illustrate this result we may take the following example. Assume that the bursting scales 
are strongly synchronized. This means that if neuron A/} ends the ith burst at a time ft, 
the neuron A4 ends the zth burst at a time t\ = tfj + £ l , where £ l ~ 0(rj) <C 0(1). Next, 
consider that within any burst in neuron A/} there are always two spikes equidistant in time. 
Let us denote r" the time at which the nth spike occurs in A/}. In neuron A4 there are two 
spikes within a burst and with a probability pk a third spike may occur [ Fig. 0(a)] . Under 
this construction, it is clear that the spiking scales are not synchronizer 
We can verify this by applying the same approach as in Refs. j^, 
threshold for the burst occurrence, the dot gray line in Fig. [7(a). Then for every burst 
we assume that the phase is increased by 2ir and between two bursts the phase increases 
linearly. So, the phase for the neuron A4 can be written as: 

«) = 2ttx + (14) 

A similar equation can be written for A/}. Note only that at a time t, the neuron A/} 

may present m bursts. So, the phase difference |A0| = \4>k{t) — <fij(t)\ is equal to \2n x 
/ t— \ / t—t m \ 

I i + i+i k i ) — 2-7T x m + +1 3 |. Now bringing the fact that \t\ — t)\ < 0(r]), we have 
V t k - l k) \ t j J 

|A0| < 4tt. (15) 

Therefore, the phase difference is bounded. On the spiking scale the situation is different; 
there is no synchronization. Doing the same procedure, we introduce a phase ip that is 
increased by 2ir between two successive spikes. Thus, the phase for the neuron A4 can be 
written as: 

A naive computation in the limit t — > oo shows that |A^| ~ pkfi. Hence, there is, of course, 
no synchronization on the spiking scale. 

Next, we construct the set Dj observing the neuron A/} every time that an event happens 
in the neuron A4- First, we fix the event to be the ending of a burst, see Fig. [7J(a). As we 
observe A/} at t\ all the points of Vj will be close to the end of the burst. So, the set Vj 
does not spread over the attractor, see the gray points in Fig. [TJb). However, the set Vj 
is also localized even if we set the event to be the occurrence of a spike. Since the spikes 
always occur within a burst, even though the spikes themselves are not synchronized, the 



Mt)=2nx (n+ m k ). (16) 
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trajectory related to the hyperpolarization period will not the visited, and therefore, the set 
T>j will be localized, see the black balls in Fig. [7(b). 
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FIG. 7: (Color online) In (a) we present the time series of the membrane potential of two neurons 
A/fc and A/j in light gray and black respectively. We show the threshold, in dashed line, for the 
burst occurrence, and in light gray dots, for the spike occurrence. While the bursting scale is 
synchronized the spiking scale is not. However, both scales can be used to construct the sets T> 
and they will be localized due to the synchronization in the bursting scale. In (b), we show the set 
T>k constructed using the spiking scale (black balls) and the set T>k constructed using the bursting 
scale (black squares). As one can see both are localized. 

VI. NEURONAL DYNAMICS 

Next, we study the appearance of PS between two spiking/bursting neurons of the 
Hindmarsh-Rose (HR) type. In such an oscillator the introduction of a phase is rather dif- 
ficult, since the neurons are non-coherent. We couple the neurons via inhibitory synapses, 
which introduces non-coherence in both time-scales. This happens because when one neu- 
ron spikes it inhibits the other neuron, which hyperpolarizes, but the neuron that has been 
inhibited still tries to spike. This competition generates even more non-coherence in both 
time-scales. Therefore, we consider this model as a proper example to illustrate our results. 
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In the 4-dimensional HR model [5|, [29( neurons are described by a set of four coupled 
differential equations: 

Xk = ay k + bx\ - cx\ - dz k + h + g S ynCl syn (x) 

Vk = e - y k + fx 2 k - gw k , (17) 

z k = fi(-z k + R(x k + H)), 

w k = u(-kw k + r(yk + l)), 

where Xk represents the membrane potential of the neuron J\f k , y k is associated 
with fast currents exchange and (zk,w k ) with slow currents dynamics, I syn (x) = 
(I syn (xi), I S yn{xi), ■ ■ ■ , I S yn(xN)) 1S the synaptic input vector and I syn (xj) is the synaptic 
current that neurons A/} (post-synaptic) injects in J\f k (pre-synaptic), and C = {%} is 
the N x N connectivity matrix where Ckj = 1 if neuron Mj is connected to neuron Nk, 
and Ckj = 0, otherwise, with j ^ k. This model has been shown to be realistic, since it 
reproduces the membrane potential of biological neurons [3^], and it is able to replace a 
biological neuron in a damaged biological network, restoring its natural functional activity 
31] . It also reproduces a series of collective behaviors observed in a living neural network 
5]. The parameters of the model are the same as in Ref. [5], but the intrinsical current Ik- 
We change Ik in order to obtain a spiking/bursting behavior and we use it as a mismatch 
parameter. First, we consider two neurons Mj and Mk- In the following, we consider the 
parameters Ik = 3.1200, Ij = 3.1205 and g syn = 0.85. 



The chemical synapses [32| are modeled by: 



[1 - Sooixi)] rS(t) = Sooixi) - S(t), 

where Xj is the post-synaptic neuron, x rev is the reversal potential for the synapse, and r is 
the time-scale governing the receptor binding. is given by: 



, tanh Xl Xth , if Xi > x th 
S OQ (x i ) = { K'^J 1 (19) 

otherwise 

The synapse parameters are x t h = —0.80, x s i ope = 1.00, x rev = —1.58. They are chosen in 
such a way to obtain an inhibitory effect in the chemical synapse. 
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FIG. 8: (Color Online) PS between two HR neurons coupled via inhibitory synapses. We analyze 
the effect of different threshold levels on the detection of PS. We depicted the attractor projection 
(x, y) in gray, and the set T> k in black, for Figs, (b) and (d). In (a) the time series of the membrane 
potentials (full lines) and the threshold Xb = —1.3 (dashed line) are depicted. In (b) the set T>k, 
construct by means of the threshold Xb, is localized; showing the presence of PS. The time series of 
the membrane potentials (full lines) and the threshold x s = 1.1 (dashed line) are depicted in (c). 
The spikes are not in PS. With our method, even for this threshold one can obtain a localized set. 
In (d), the set Uj- construct by using the threshold x s = 1.1 is localized. 

To construct the sets T>, we define the event occurrence. We shall analyze two situations: 
when the event is defined in the bursting scale, and when the event is denned in the spiking 
scale. Firstly, we define the iih event to be the ith crossing of the membrane potential of 
the neuron Mj t k with the threshold Xb = —1.3 in an upwards direction. We denote the time 
events by ft k . Note that this threshold assigns to the times t % - k the beginning of the zth 
burst of Afj,k- Fig. [8](a) shows the time series of the membrane potential of the neurons Mj,k- 
The threshold Xb = —1.3 is depicted with the dashed line, and it is chosen in such a way 
that it does not define a proper Poincare section, which means that not all the bursts cross 
it, see FigE^a). Actually many bursts are missed. Thus, the approach to extract the phase 
considering the increasing of 2tt between two bursts, misleads the statement of PS. That 
is so, because in this approach the phase is threshold dependent. Therefore, by using Eq. 
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(|T4|) . we get that PS does not exist, which is crucially wrong (note that, with the increasing 
of the threshold value PS would appear). However, our approach, which is not threshold 
dependent, overcomes these difficulties. Indeed, localized sets exist even for this threshold 
[Fig Kb)]. 

Conversely, if we increase the threshold level in such a way that it takes into account 
the spike occurrence, e.g. a threshold at x s = 1.1, the dashed line in Fig. [8](c), the former 
approach, as in Eq. ffT6]) . completely fails to state PS, due to the fact that the spikes are 
not in PS. Furthermore, the spikes are highly non-coherent. The competition between the 
two neurons generates a damping in the spikes in the beginning of the burst, followed by an 
increasing and then decreasing in the spike frequency [Fig. |Sfc)]. Again, since the threshold 
x s = 1.1 defines a typical event, the observation of A/},* at times (t\ j)ieN provides localized 
sets V [ Fig. MM- 

VII. EXCITATORY NEURAL NETWORKS 

The ideas introduced herein are also useful to analyze the onset of synchronization in 
networks. We consider a network of 16 non-identical HR neurons, regarded as Mi where 
% G [1, . . . , 16], connected via excitatory chemical synapses. The mismatch parameter is 
the intrinsic current l{. Since the meaningful parameter is Ii = 3.12, for which the HR 
neuron best mimics biological neurons, we introduce mismatches around this value for all 
the neurons within the network. Thus, given a random number 77, uniformly distributed 
within the interval [—0.05,0.05], we set Jj = 3.12 + rji. The excitatory synapses are modeled 
by Eqs. (18) and ([19]) . To obtain the excitatory effect we change the value of x rev . If 
pre-synaptic neuron always injects a positive current in the post-synaptic 
one. Since the maximum spike amplitude is around 1.9, we set x rev = 2.0. 

Our network is a homogeneous random network, i.e. all neurons receive the same number 
k of connections, namely k = 4, see Fig. [9](a). We constrain g syn [see Eq. (17)] to be equal 
to all neurons. We identify the amount of phase synchronous neurons by analyzing whether 
the sets T>i are localized. 

The onset of PS in the whole network takes place at g* yn ~ 0.47; so all neurons become 
phase synchronized. As the synapse strength crosses another threshold, g syn ~ 0.525, the 
neurons undergo a transition to the rest state no longer presenting an oscillatory behavior. 



18 



(a) (b) 

FIG. 9: Networks generated randomly. In (a) n=16 and k=4, while in (b) n=9 and k=3. 

Clusters of PS appear for g syn <C g* yn . In fact, right at g syn ~ 0.04, some PS clusters appear [ 
Fig. HOTa)]. Again, the clusters are identified by analyzing the localized sets. These clusters 
seem to be robust under small perturbations. 

Clusters of PS inside the network may offer a suitable environment for information ex- 
changing. Each one can be regarded as a channel of communication, since they possess 
different frequencies, each channel of communication operates in different bandwidths. To 
see the bandwidths in the network, we analyze the variance in the average bursting time of 
the neurons. Since only the burst scale is synchronized, we are just interested in the average 
bursting time, which can be straightforwardly estimated with a fast Fourier transformation 



FFT 33|. So, given the neuron A/j, we label its bursting average time by (Tj). Then, we 
compute the variance of the average time on the ensemble of neurons. For this, we first 
introduce the average time of the whole network, which is given by: 



1 n 

c = -£&■>■ ( 2 o) 



n 

3=1 



Thus, the variance of the average time on the ensemble of neurons is readily written as: 

1 71 

- = -E(^)-o 2 - (2i) 



n . 



So, a indicates how diverse are the bandwidths. As one can see in Fig [TOTb). when the 
first clusters appear for g syn w 0.04, we have a ~ indicating that the whole network is 
working almost with the same frequency. A further increasing of g syn causes the destruction 
of these clusters and an increasing of a. However, even in the regimes of high a with 
Qsyn £ [0.27, 0.34], there is the formation of clusters. 
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FIG. 10: The appearance of PS clusters within the network. In (a) we show the number of clusters 
as a function of the synaptic strength. In (b) we plot the normalized a*, where a* = cr/112.8. 
Note that when the clusters are formed a becomes small, but bounded away from zero, which 
means that the neurons within the network undergo a transition where they have almost the same 
frequency. 

This scenario of cluster formation is neither restricted to this HR model nor to the 
synapse model. It can also be found in square-wave and parabolic bursters, and it is in 
general achieved quite before the onset of complete synchronization. For example, we use a 
more simplified HR model given by: Xj = ax 2 — x^ — yj — Zj — g syn (xj — x s )Cl syn ('x),y = 
(a + a)x 2 — y, z = fi(bx + c—z), with the parameters: a = 2.8, a = 1.6, c = 5, b — 9, // = 0.001; 
C being the connectivity matrix and I SJ)n (x) = (I syn (xi), . . . , I S yn(%N)) a fast threshold 
modulation as synaptic input given by 

I syn (xj) = 1/[1 + exp{-(3{ Xj - 6)}], (22) 

with (3 = 10 and O = —0.25. As before, g syn is the synaptic strength and the reversal 
potential x s > Xj(t) in order to have an excitatory synapse. For a homogeneous random 
network of 9 identical HR neurons, with k = 3 [ Fig. [9]^b)], the theory developed in Ref. 



341 ] predicts the onset of complete synchronization at g syn ~ 0.425, while we found that PS 



in the whole network is already achieved at g* yn ~ 0.36. Clusters of PS, however, appear 
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for a much smaller value of the coupling strength, actually at g syn ~ 0.03. Next, we apply 
the same procedure as before and we compute the variance of the average bursting time on 
the ensemble of neurons within the network. The result a x g syn is depicted in Fig. [TTJ the 
inset numbers indicate the amount of clusters. 




FIG. 11: The average bursting time on the ensemble of neurons. We plot the normalized a*, 
where a* = a/12. The inset numbers show the amount of clusters for a given parameter g syn 



As we have pointed out, such clusters are rather suitable for communication exchanging 
mainly for two reasons: (z) They have different frequencies, therefore, each cluster may be 
used to transmit information in a particular bandwidth, which may provide a multiplex- 
ing processing of information, (ii) The clusters of phase synchronous neurons provide a 
multichannel communication, that is, one can integrate a large number of neurons (chaotic 
oscillators) into a single communication system, and information can arrive simultaneously 
at different places of the network. This scenario may have technological applications, e.g. 
in digital communication [37], |38[ , and it may also guide us towards a better understanding 



of information processing in real neural networks jl5|, lly, l39j ]. 



VIII. DETECTION OF PS FOR HIGHER DIMENSIONAL SYSTEM 

It is easy to say whether the set T> is localized in a two dimensional plane; this could 
be done for example by visual inspection. In mult i- dimensional system it might not be 
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obvious whether the set T> is localized. This is mainly due to the fact that in a projection 
of a higher dimensional system onto a low dimensional space, the set T> might fulfill the 
projected attractor. Therefore, the analysis of the localization might have to be realized in 
the full attractor of the subsystem. 

The analysis is also relatively easy if we bring about a property of the conditional ob- 
servation. Whenever there is PS, the conditional observation, given by F fe (x°), is not 
topologically transitive [3^] in the attractor of £j, i.e. T> is localized. The conditional ob- 

4-i I I 

servations F J fe (x°) are topologically transitive in the attractor Aj of £j [36J if for any two 
open sets B, C C A, 



3t n> j F h (B)nC^ 0. (23) 



To check whether T>j is localized, we do the following. If there is PS, for e T>j it exists 
infinitely many Xj G Aj such that 

Yj n Bfa) = 0, (24) 

where Bt(xj) is an open ball of radius I centered at the point Xj, and i is small. We may 
vary yj and Xj to analyze whether it is possible to fulfill Eq. ff24|) . Whenever this is possible, 
it means that the set Vj does not spreads over the attractor of and therefore, there is 
PS. 



For analysis of PS basing on experimental data [111, ll8|, |20| where the relevant dynamical 
variables can be measured, so that the phase space is recovered, our approach can be used 
straightforwardly. If one just has access to a bivariate time series, one first has to reconstruct 
the attractors, and then proceed the PS detection by our approach 

IX. INFORMATION TRANSMISSION 

In this section, we analyze the relationship between the sets T> and the capacity of infor- 
mation transmission between chaotic oscillators. In order to proceed such an analysis, we 
may assume that the oscillators are identical or nearly identical. Such that the synchronized 
Tajectories are not far from the synchronization manifold, i.e. the subspace where Xj = X& 
3j. Next, for a sake of simplicity we consider only oscillators whose trajectory possess a 
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proper rotation and are coherent [22J, |38|, e.g. the standard Rossler oscillator. However, the 
ideas herein can be extended to other oscillators as well. 

The amount of information that two systems Ej and E^ can exchange is given by the 
mutual information J(E 3 -, E^) j^oj]: 

I(X j ,E k ) = H(X j )-H{E j \X k ), (25) 

where H(Ej) is the entropy of the oscillator Ej and H(Ej\E k ) is the conditional entropy 
between Ej and E^ , which measures the ambiguity of the received signal, roughly speaking 
the errors in the transmission. 

As pointed out in Ref. 14j the mutual information can be also estimated through the 
conditional exponents associated to the synchronization manifold. The mutual information 
is given by: 

7(E;,E fc ) = ^A+-J>+ (26) 

where Ay" are the positive conditional Lyapunov exponents associated to the synchronization 
manifold, the information produced by the synchronous trajectories, and \~\_ are the positive 
conditional Lyapunov exponents transversal to the synchronization manifold, related with 
the errors in the information transmission. In PS A J can be small, which means that one 
can exchange information with a low probability of errors. So, PS creates a channel for 
reliable information exchanging 14j. In general, we expect ^ Ajj~ < where A + are 



the positive Lyapunov exponents. Thus J(Ej,Efc) < \ + — ^ A j\ In order to estimate 
an upper bound for /(Ej, E^), we need to estimate Xj_, what can be done directly from the 
localized sets. 

The conditional transversal exponent can be estimated from the localized sets by a simple 
geometric analysis. At the time the oscillator Ej reaches the Poincare plane at x* while the 
oscillator E^ is at x.\ = x& (£*•). The initial distance between the trajectories is Ax^ = x*— x^.. 
This distance evolves until the time t\ when the oscillator E& reaches the Poincare plane 
at x* k , while the trajectory of Ej is at xj = jc.j(t l k ). The new distance is Ax.jk(t l k — ft) = 
x* k — ~x.j(t k ). Therefore, we have: 

Ax jfc = Ax ife e A il^-*il (27) 
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So, the local transversal exponent is given by: 

1 
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Xj x A 



(28) 



where we use the convention x logO = 0. Of course, we only estimate the conditional 
exponent close to the Poincare plane. Hence, if we change the Poincare plane the conditional 
exponent may also change, i.e. there are some events that carry more information than 
others. 



A. Example with Rossler Oscillators 

We illustrate this approach for two coupled Rossler oscillators. We set the parameters 
to a = 0.15, b = 0.2, c = 10, a s = 1, and Aa k = 0.0002. As shown in Ref. Q at 
e w 0.05, the two oscillator undergo a transition to PS. In particular, for e = 0.06 we 
have Yl ^± ~ 0.06. We estimate a ^ this situation by means of Eq. (128]) . We set 

the Poincare section at yj^ = 0, and compute Aj_ for 65, 000 cycles, i.e. 65,000 crossing 
of the trajectory with y = and y < 0. We get Aj_ ~ 0.048. Note that we are not 
computing J^A^, but rather, the maximum A^, namely Aj\ Therefore, it is natural to 
expect A^ to be smaller than Y2 However, the upper bound to the information exchange 
can be estimated by I(S,R) < Y A + — Aj, that is, the maximum amount of information 
that can flow through the coupled oscillators if we encode the trajectory using the Poincare 



plane y = 37|]. Furthermore, it seems that when the level of synchronization is large, 
the estimation of A J, by means of Eq. (128|) . might become problematic, due to strong 
fluctuations in \t) - f fc | _1 in[\(x* - xj; )/(x£ - x*)|]. 



X. CONCLUSIONS 



We have proposed an extension of the stroboscopic map, as a general way to detect PS 
in coupled oscillators. The idea consists in constraining the observation of the trajectory 
of an oscillator at these times in which typical events occur in the other oscillator. This 
approach provides an efficient and easy way of detecting PS, without having to explicitly 
calculate the phase. We have shown that if PS is present, the maps of the attractor appear 
as a localized set in the phase-space. This has been illustrated in coherent oscillators, the 
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coupled Rosslers, as well as in non-coherent oscillators, spiking/bursting neurons of HR type 
coupled with chemical synapses. As we have shown in neural networks, the appearance of 
clusters of PS is rather common, which may be relevant for communication mainly due to two 
aspects: (i) The clusters provide multiplexing information processing, namely each cluster 
may be used to transmit information within a bandwidth, (ii) They provide a multichannel 
communication, that is, a large number of neurons is integrated into a single communication 
system. Moreover, we have analyzed the relation between the information exchanging and 
the localized sets. We have roughly estimated the errors in the information transmission 
from the localized sets. 
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APPENDIX A: PROOF OF THEOREM 1 

In this appendix we prove the theorem 1. It is instructive to give a sketch of the proof, 
in order to have a better understanding of the result. We split the demonstration into the 
following four steps: (i) We show that the increasing of 2ir in the phase 0^ defines a smooth 
section Yj^ on Ej^, which does not intersect itself, (ii) We show that observing the oscillator 
Sj whenever oscillators crosses Yk gives place to a localized set Vj. (Hi) Further, we 
show that the observation of Sj whenever S fc crosses a piece Pr k of the section Y k also gives 
place to a localized set. (iv) Using these results we show that, actually, the localized sets 
can be constructed using any typical event. To show this, we only note that given a typical 
event with positive measure, we can choose Pr k to be close to the event occurrence, implying 
that shortly before or shortly after of every event occurrence, a crossing of the trajectory 
with Pr fe will happen. Thus, if we observe Ej whenever the event occurs in we will have 
a set that a close the T> j: and therefore, localized. Next, we formalize the heuristic ideas. 
Let us introduce i = j, k. 

Proposition 2 The increasing of2n in <j>i(t) generates a smooth section Tj in the attractor 
of Hi, which does not intersect itself. 
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Proof: Firstly, let us introduce the times (t™ - ) such that 4>i{T™) = m x 2tt. Then, let Tj 
be the set of points such that given the initial point we have the section: 

r, = {u meN xn xr = Ff (xf)} (Al) 

Thus, we construct a section Tj. Ti is smooth since both 0, and F* are smooth. Indeed, 
given two points x°, xj G r i; with d(x°, x*) < e, there is a r > 1 such that (x°), F^ (x ? - ) G 
Tj, and 

d(F?(x°),Ff(x l 1 ))<5. (A2) 
Furthermore, we can construct a continuous section Tj, by conveniently choosing points 



xf. The fact that Tj does not intersect itself comes from the uniqueness of F* [4l||, and from 
the fact that the <j>i(t) > 0, which implies that the phase is an one-to-one function with the 
trajectory Note that, obviously, this section depends on the initial conditions. □ 

Lemma 1 The observation of the oscillators whenever the trajectory of crosses the 
section T k gives place to a localized set Vj if, and only if, there is PS. 



Proof: Let II, be the Poincare map associated to the section Tj, such that given a point 
x™ G Tj, so x" +1 = ilj(x^) = 3 ( x j)? where At™ = r™ — t™ -1 . From now on, we use a 
rescaled time t' = t/(Tj), with (Tj) = linij_ >0O Tj/z. For a slight abuse of notation we omit 
the /. There are numbers Ki such that \t- — i(Ti)\ < Ki, where, by time reparametrization, 
Kj<l. If both oscillators are in PS, then (T&) = (Tj), and so: 

K-r-\<k, (A3) 

with k < Kfc + Kj <C 1 . Now, we analyze one typical oscillation, using the basic concept 
of recurrence. Given the following starting points x° G Tfe and x° G Tj, we evolve both until 
x° returns to Tj. Let us introduce 

Ar n = At™ - At™. (A4) 

Which gives: 

F^(x°) = n,(x°)=x)Gr,, (A5) 

Analogously, 
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Ff.'(x») = Ff^'K) 

= Ff'oFf'tx"). 

Bringing the fact that F k Tk (x°) = IT fc (x°) = x\, we have: 

Pf^(xS) = Ff' a (xi). (A6) 

Now, by using the fact that |Ar*| < k, we can write: 

Ff 1 (xi)«xi+G(x^ + 0(K 2 ). (A7) 

So, given a point x/. G evaluated by the time when the trajectory of Ej returns to 

the section Tj, the point x& returns near the section T^, and vice-versa. Therefore, it is 

localized. For a general case, we have to show that a point, in the section T^, evolved by the 

flow for an arbitrary number iV of events in the oscillator Ej, still remains close to T^, in 

other words, it is still localized. This is straightforward, since | YliLo ^ T% \ = \ T k ~ T j^\ < ^- 

So, we demonstrated that the PS regime implies the localization of the set T> k . 

Now, we show that the localization of the set T> k implies PS. Supposing that we have a 

localized set V k , so, Eq. (1A3I) is valid, by the above arguments. Therefore, we just have to 

show that Eq. (IA3I) implies PS. With effect, we have \(fij(t) — (fik(t)\ — I f ^jdt — f Q fl k dt\ 

which is equal to | J J Qjdt — J J Q k dt + J n Qjdt — J n Q k dt\. This may be written as 

| Jq j Qjdt — J Tfe Q k dt — J T ^ n k dt + f l „ Qjdt — /*„ Q k dt\. Next, noting that 0j(r") = 2n x n, 

T j T j T j 

we get: 

- Mt)\ < M\t? - r£| + 2AM, (A8) 

where A = max \t? - Therefore, if the t ime event difference — f%\ is bounded it 

implies the boundedness in the phase. Thus, we conclude our result. □ 

Proposition 3 Let (tJ**)^^ ^ e the times at which the trajectory Ej crosses a piece Py 3 of 
Tj. If there is PS, then the observation of the trajectory ofE k at times (rj 1 *)^^ gives place 
of a localized set. 

Proof: Note that the observation of the trajectory of E^ at times (tj)^ gives place to a 
set V k , while the observations at times (r^)^^ give place to a subset V k of V k . Therefore, 
whenever T> k is localized, it implies the localized of T> k . □ 
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Now, we are ready to prove the theorem 1. 

Proof: Let the event be the entrance in an e-ball, such that the event occurrence produces 
the time series tj, in Ej. There is, at least, one intersection of this ball with the section Tj. 
Since Tj depends on the initial conditions, we can choose an initial condition right at the 
£-ball event. Next, we choose P? such that it is completely covered by the £-ball. Since the 
measure of the £-ball is small, s < 1, the time difference between crossings of the trajectory 
with Pp and the £-ball is small, thus, there is a number rj < 1 such that: 

t) - 7f < 0(?7) • (A9) 

Therefore, if we observe the trajectory of E^ at times (£})ieN> we have a localized set in E/-. 
Thus, we conclude our result: The observation of the trajectory of E^j. whenever typical 
events in E^j occurs generates localized sets "Dj^ if; an d only if, there is PS. □ 
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